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SUMMARY 

The  author  has  previously  used  the  stochastic  model  of 
radiative  transfer  to  obtain  the  diffusion  matrix  consisting 
of  the  reflectance  and  the  transmittance  operators  In  a  finite 
plane— parallel  atmosphere  of  arbitrary  stratification.  When 
this  method  Is  applied  to  the  distribution  of  emission  sources 
within  the  medium,  the  emergent  Intensities  can  be  computed 
without  actually  solving  the  equation  of  transfer. 

In  the  present  paper,  using  auxiliary  equations  in 
conjunction  with  the  Milne  integral  equations,  the  author 
derives  the  diffusion  matrix  along  with  the  extension 
concerning  the  Neumann  solution  aB  given  by  Busbrldge.  In  the 
case  of  diffuse  reflection  and  transmission  of  parallel  rays, 
the  solutions  are  expressed  in  terms  of  a  pair  of  scattering 
and  transmission  functions  for  each  of  the  two  boundaries  of 

# 

the  atmosphere.  Then  these  global  functions  are  given  by  X 
and  Y*  functions  that  are  equal  to  those  previously  found  by 
Bellman  and  Kalaba.  Whereas  the  diffusion  matrix  formally  has 
a  somewhat  similar  appearance  to  a  map  yielded  by  Preisendorfer, 
the  mathematical  development  is  different. 

If  the  optical  properties  of  the  medium  are  constant 
throughout  the  atmosphere,  the  reflectance  and  transmittance 
operators  obtained  here  reduce  to  those  given  by  Sobolev. 
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*  ON  THE  DIFFUSION  MATRIX  OF  RADIATIVE  TRANSFER 
I.  INTRODUCTION 

In  1942  Ambarzumlan  [l]  published  a  mathematical  method 
for  solving  the  homogeneous  first  integral  equation  of  Milne. 
Basing  his  work  on  the  idea  of  linear  aggregation,  he  reduced 
the  solution  of  the  Milne  equation  to  that  of  its  auxiliary 
equation,  which  in  turn  leads  directly  to  thu  H  equation  (see 
Kou'rganoff  [22]).  in  a  second  paper,  Ambarzumlan  [2]  treated 
the  same  problem  from  a  physical  viewpoint.  By  means  of  the 
principle,  of  invariance,  he  obtained  a  nonlinear  integral 
equation  in  the  H  function  directly  through  physical  analysis 
of  the  diffuse  reflection  and  transmission  of  light,  without 
actually  solving  the  homogeneous  Milne  equation  for  the  source 
function.  Ingeniously  extending  the  invariance  method, 
Chandrasekhar  [15 ]  formulated  a  complete  set  of  principles  of 
invariance  in  a  finite  homogeneous  atmosphere  and  applied  it 
to  various  transfer  problems.  The  requisite  integral  equations 
for  the  scattering  and  transmission  functions  were  derived  from 
the  above  set  of  invariance  principles  in  connection  with  the 
equation  of  transfer. 

By  means  of  the  Ambarzumlan  technique,  the  solution  of  the 
nonhomogeneous  Milne  equation  can  be  expressed  in  terms  of  the 
H  function  fox*  certain  special  forms  of  the  distribution  of 
emission  sources  acting  in  the  medium.  The  technique  plays  an 
important  role  in  the  solutions  of  transfer  problems  from  the 
practical  viewpoint.  Furthermore,  as  Busbridge  [?„3 ]  emphasized. 
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the  technique  becomes  more  powerful  when  it  is  connected  with 
the  theory  of  the  N  solution.  The  auxiliary  equation  for  any 
given  problem  reduceB  to  the  first  integral  equation  of  Milne 
for  the  diffuse  reflection  and  transmission  of  a  parallel  beam 
of  radiation  by  a  finite  layer.  Prom  the  mathematical  point  of 
view,  the  Ambarzumian  technique  was  Ingeniously  extended  by 
Busbrldge  [12,13]  and  was  applied  to  various  transfer  problems 
of  aBtrophysical  Interest  by  Busbrldge  [12,13 j  and  Ueno  [29 J. 

Ambarzumian 1 s  physical. method  based  on  the  principle  of 
Invariance,  has  been  applied  by  Busbrldge  [11,14]  and  Stibbs 
[14,25]  to  problems  of  line  formation  in  the  Mi'lne-Eddington 
model  with  coherent  and  noncoherent  scattering. 

Recently,  developing  the  idea  of  the  invariance  principle 
of  Ambarzumian  [2],  Bellman  and  Kalaba  [3]  stated  the 
principle  of  invariant  imbedding.  The  functional  relationships 
among  the  members  of  the  class  are  found  by  imbedding  the 
original  process  within  a  family  of  processes  of  similar  nature 
and  obtaining  an  Invariant  process.  Whereas  the  classical 
approaches  reduce  problems  to  the  solutions  of  Bystems  of 
linear  equations*  the  invariant-imbedding  technique  reduces 
problems  to  the  iteration  of  nonlinear  transformations.  The 
principle  of  invariant  imbedding  led  not  only  to  new 
analytical  functions  of  radiative  transfer  (Bellman  and  Kalaba 
[2,3]),  neutron  diffusion  (Bellman,  Kalaba,  and  Wing  [8,10]), 
random  walk  and  scattering  (Bellman  and  Kalaba  [6]),  adaptive 
processes  and  random  transmission  (Bellman  and  Kalaba  [7]),  and 
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problems  of  Stefan  type  (Bellman  and  Kalaba  [9l),  but  also  to 
those  of  wave  propagation  (Bellman  and  Kalaba  (5))* 

By  means  of  the  Chandrasekhar-Wick  method,  and  with  the 
aid  of  Chandrasekhar's  extension *of  the  physical  method  based 
on  the  invariance  principle,  the  following  kinds  of  transfer 
problems  have  been  solved s  line  formation  in  planetary  atmos¬ 
phere  (Chandrasekhar  [15 1 ) *  line  formation  in  the  Schuster 
model  (Chandrasekhar  [15 3) *  emitting  atmosphere  (Horak  [l8], 
Horak  and  Lundquist  [19]),  and  molecular  absorbing  atmosphere 
(King  [21]).  . 

Allowing  for  the  map  consisting  of  complete  reflectance 
and  complete  transmittance  operators  in  inhomogeneous  one- 
parameter  carrier  space,  Preisendorfer  [23  ]  obtained  the 
invariant  imbedding  relation  for  radiative  transfer  and 
neutron  transport  contexts.  Furthermore,  the  functional 
relations  for  reflectance  and  transmittance  operators  were 
derived  in  a  manner  similar  to' that  used  by  Chandrasekhar  for 
the  homogeneous  case  [24].  The  polarity  of  the  two  operators 
is  elucidated  because  of  the  inhomogeneous  optical  properties 
of  the  medium. 

In  recent  years,  introducing  the  probability  concept  into 
the  theory  of  radiative  theory,  Sobolev  127)  has  treated  many 
subjects  from  the  statistical  point  of  view,  including  various 
kinds  of  transfer  problems:  pure  scattering  [27),  line 
formation  with  coherent  and  noncoherent  scattering  in  a  semi- 
infinite  homogeneous  medium  [27],  diffuse  reflection  in  a 
semi— infinite  inhomogeneous  medium. [26],  diffuse  reflection  in 
a  finite  homogeneous  medium  [28],  and  others  [27  ]. 
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Recently,  assuming  that  multiple  scattering  of  a  photon 
as  the  carrier  of  radiant  energy  is  a  random  process  of 
Markovian  type,  Ueno  used  a  stochastic  approach  in  the  study 
of  various  problems  of  radiative  transferj  Milne's  problem 
[JO],  line  formation  in  semi- infinite  atmospheres  with  coherent 
and  noncoherent  scattering  [31,32],  diffuse  reflection  and 
transmission  in  a  finite  homogeneous  and  inhomogeneous  layer 
[33*3^1*  and  a  Markovian  property  of  radiative  transfer  [35]. 
Through  the  probabilistic  technique,  the  integro-differential 
equation  for  the  emission  probability  distribution  can  be  ' 
derived  from  the  Chapman— Kolmogoroff  equation,  and  the  integral 
transform  of  the  probabilistic  equation- can  be  used  to  obtain 
the  angular  distributions  of  the  emergent  radiations  without 
actually  solving  the  equation  of  transfer.  •  It  is  of  interest 
to  mention  that  the  derivation  of  the  S  and  T  functions 
from  the  probabilistic  equation  by  the  stochastic  approach  is 
similar  in  part  to  that  used  recently  by  Feller  [16]  on 
boundaries  and  lateral  conditions  for  the  Kolmogoroff 
differential  equation. 

In  a  preceding  paper  [37],  based  on  a  stochastic  model  of 
radiative  transfer,  the  author  obtained  the  diffusion  matrix 
consisting  of  reflectance  and  transmittance  operators.  In  the 
present  paper,  however,  we  derive  the  diffusion  matrix  from  the 
auxiliary  Equations  that  play  a  major  role  in  the  Ainbarzumiari 
technique  in  cqnnection  with  the  Milne  equations  (cf.  Hopf  [17], 
Bu3bridge  [13]).  In  subsequent  papers,  the  diffusion  matrix 
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will  be  applied  to  the  various  transfer  problems  of  current 
interest  in  a  finite  inhomogeneous  medium. 

Whereas ,  in  the  context  of  the  operational  form,  the  map 
due  to  Preisendorfer  [23]  seems  to  be  somewhat  similar  to  the 
diffusion  matrix,  the  mathematical  procedure  is  different.  In 
contrast  to  Preisendorfer ' s  approach,  which  can  be  patterned 
after  the  p  rocedure  used  by  Chandrasekhar  for  a  finite  homo— 
geneous  atmosphere,  our  technique  is  based  on  the  uBe  of 
auxiliary  equations  for  the  source  functions. 

Finally,  it  should  be  mentioned  that,  in  addition  to  the 
methods  stated  above,  the  Laplace— transform  method  (see  Hopf 
[17],  Kourganoff  [22],  Huang  [22  ], . and  BuBbridge  [13])  is  also 
used  in  transfer  problems . 

II.  THE  EQUATION  OF  TRANSFER 

Let  Iv(tv,h)  be  the  specific  intensity  of  v  radiation 
at  the  optical  depth  Ty  in  the  direction  cos- V.  In  con¬ 
sidering  the  transfer  problem  of  radiation  in  a  plane-parallel 
inhomogeneous  atmosphere  of  finite  optical  thickness  ,, 
for  simplicity  we  shall  restrict  our  discussion  to  coherent 
and  isotropic  scattering.  Then,  for  convenience,  we  suppressed 
the  subscript  v  of  the  various  quantities. 

The  equation  of  transfer  appropriate  to  the  present  case 
is  written  in  the  form 

(2.1)  |l  .  !(•,,„)  -  j(,), 


where  the  source  function  J(t)  is 
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(2.2)  J(t)  -  £  w(t)  J  l(T,H')d|A»  +  B1(t). 

-1 

In  equation  (2. 2),  w(t)  ia  the  albedo  for  single  scattering, 
and  B^(t)  is  the  distribution  of  emission  sources  within 
the  atmosphere. 

Following  the  notation  of  Chandrasekhar  in  [15 )»  we  shall 
denote  the  intensity  of  radiation  at  the  level  t  directed 
toward  the  surface  t  -  0  by  I(t,+  h),  0  <  H  £  1,  and  that 
directed  toward- the  surface  t  -  by  I(t,—  |i),  0  <  |a  £  1. 

Let  radiation  of  intensity  1(0,-  n),  0  <  |a  £  1,  be  Incident 

* 

on  the  surface  t  •>  0.  in  the  direction  —  n,  and  let  radiation 
of  intensity  I(t1#+  h),  0  <  H  £  1,  be  incident  similarly  on 

the.  surface  t  -  t1  in  the  direction  +  |x. 

Then,  allowing  for  the  formal  solution  of  equation  (2.1) 
subject  to  the  boundary  conditions  given  above,  namely’, 

(2.5)  I(t,+  |x)  a  i(t1,+  n)exp( - - ) 

+  /Tl  J(t)exp<- 

T 

(2.4)  I(t,-  H)  -  1(0,-  n)e*p(-  1) 

+  /  j(t)exp(- 
0 

from  (2.1)  we  get  the  first  integral  equation  of  Milne  ' 
governing  J(t)j 
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(2-5) 

[1  -  w(T)A]T(j(t))  =  B(t), 

where  1 

is  the  identity  operator,  H  is  the  truncated  Hopf 

operator 

• 

(2.6) 

V'(»M  -  1/ 

Tl  f(t)?:1(|t  -  t | ) dt , 

0 

and  B(t) 

is 

• 

• 

(2.7) 

B(t)  -  Bx(t)  + 

I  w(t)  J  1(0,- n')exp(- jjr)dn' 

* 

0 

+  §  w(t)  J 

^  I(i1,+  H')exp(-  \,-I)d,.. 

0 

-  ♦  • 

In  equation  (2.5),  E1(t)  is  the  first  exponential  integral: 


(2.8)  Ei(t)  =  J  1  exp(-  I)f.. 

0 

While  in  the  homogeneous  case  B(t)  is  a  known  function 
in  any  particular  problem,  the  quantity  B^t)  in  the  case  of 
inhomogeneous  distribution  of  emission  sources  is  given  by 

(2.9)  B^(t)  '=>  J  1  ^  J'1  B,(t,+  u)exp(-  ^-=-^)dt 

0  L  T 

+  J  B1(t,—  M-)exp(—  ^-~-)dt  . 

0 

III.  THE  AUXILIARY  EQUATIONS 

In  the  theory  of  radiative  transfer  based  on  the 
Ambarzumian  technique,  the  auxiliary  equation  plays  an 
important  role . 
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In  thb  case  of  a  finite  inhomogeneous  atmosphere,  the 
auxiliary  equations  are  expressed  in  the  forms 


(J.l)  [1  -  w(T)X]<rJp(njt,T1)}  -  w(*r)exp(-  J), 

T  T 

(3.2)  [\  -  w(T)I]<r(p*(^jt,T1)}  -  w(x)exp(-  - --jj--  ) . 


Physically,  the  function  p(m>JT,t^)  corresponds  to  the 
source  function  due  to  a  parallel  beam  of  radiation  incident 

*,  V 

on  the  surface  t  =  0  in  the  direction  —  M-,  and  p 
corresponds  to  that  due  to  similar  parallel  rays  on  the  sur¬ 
face  t  -  Tj^  in  the  direction  4  M-.  Prom  the  probabilistic 
viewpoint,  the  function  p(M';*r,T1)  represents  the  probability 
that  a  photon  absorbed  at  the  level  t  will  reappear  in  the 
direction  4  M-  as  radiation  emerging  from  the  surface  t  -  0, 
and  p*(M-jt,t1)  is  the  probability  that  a  photon  absorbed  at 
the  level  t  will  be  reemitted  in  the  direction  -  |i  as 
radiation  escaping  from  the  surface  t  ■  Mathematically, 
the  probability  p(njT,T^)dp.  1b  the  probability  of  finding  |A 
in  the  range  (^,p.  4  M-)  at  the  level  t. 


On  differentiating  the  auxiliary  equations  with  respect 

to  T^,  after  some  argument  we  have 

ap(H;T,T,)  1  r  1  *  %du, 

(3.3)  — - £  J  P  (n.;T,Tl)^T, 

1  0 

5p*(HiT,T.)  1  •* 

- 57- - p:  P  <^5T,Tl) 

+  I  J  p*(n»;T,T1)^{-. 


(5.4) 
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Equations  (3-3)  and  (5-4)  can  also  be  obtained  by  means  of  the 
probabilistic  method  [34]  and  the  invariant  imbedding  technique 
[36].  Furthermore,  the  solutions  are  found  to^depend  on  four 

functions  X(|a,t1).,  Y(d,'i1)J  *  and  Y  (^»Ti)»  *"^h 

are  connected  to  each  other  by  the  principle  of  reciprocity 
(dee  [?4]).  Mathematically  speaking,  the  source  functions 
p(M.;t,t1)  and  pVjt,^)  are  .the  N  solutions  of  the 
auxiliary  equations  (j.l)  and  (3-2)  (see  Busbridge  U3))- 


IV.  diffusion  matrix 

With  the  aid  of  equations  (2.3)  and  (2.4),  the  emergent 
intensities  from  the  atmosphere  l(0,+  M-)  and  l(Ti»“  I1)  ai 
respectively  provided  by 

(4.1)  1(0, +  n)  «  I(tx,+  H)exp(-  -jp)  +  l*(0,+  t0» 

^1  # .  » 

(4.2)  *)  -  I(°»“  ^)exP(“  +  1  *T1*“ 

where 


(4.3) 

(4.4) 


I#(0,+  H)  -  /Tl  J( t)exp(— 

0 

•  p  1  T  m  *“  t 

I#(t1,-  H)  -  J  1  J(t)exp( - jj - )-jT 

0 

./Tl  J(t1  -  t)exp(- 
0 


Following  Busbridge  [13 we  call  equations  (4.3)  and  (4.4) 
the  Tj  transforms  of  J(t). 
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When 

(4.5)  B(t)  -  aBX(T)  +  bB?(a1  -  t),- 
and 

(4.6)  [1  -  w(T)I]T(J1(t))  -  B1(t), 

(4.7)  [1  -  w(t)X]T(J2(t))  -  B2(T]L  -  t), 

then  the  solution. of  equation  (2.5)  1b  written  In  the  form 

(4.8)  J(t)  -  aJ1(t)  +  bJ2(x) . 

Furthermore,  the  emergent  Intensities  I*(0,+  p.)  and 
I  (t^,—  u)  are  given  by 

(4.9)  I*(0,+  ,)-a/Tl  J^(t)exp(—  £)££ 

0 

+  b/Tl  J2(t)exp(—  ±)f, 

0 

(4.10)  I*(Tr-n)  -  a/’1  'j’(t)exp(- 

0 

+  tJexpt-V1)"- 
0 

On  multiplying  equation  (4.6)  by  p(M-;t,t^)/w(t)m-  and 
integrating  with  respect  to  t  over  (0,t^),  we  get 

(4.11)  f1  Ip(«;t.-'1)4[t}  " 

0 
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On  the  other  hand.  If  we  multiply  equation  (3.1)  by 
J1(t)/w(t)m-  and  integrate  with  respect  to  t  over  (0,t^), 
we  get 

(4.12)  p{(MT,Ti)  _  J1(T)XTtp(M.,t,T1)J 

o 

-  f1  J1(t)cxP(— ■  !>£. 

'o 


Allowing  for  the  symmetrical  property  of  the  3C  operator 
(cf.  Busbridge  [13]),  we  have 

f  1  fi(t)fit{f2(t'))dt  »  rTl  f2(t)Xt{f1(t')}dt, 

0  0 


and  using  equations  (4.11)  and  (4.12),  we  obtain 

(4.13)  p  J1(T)exp(-  1)31  .  p  pOuT.tjjSifl}  31. 


Similarly,  using  equations  (3.1),  (3.2),  (4.6),  and.  (4.7), 
we  get 


(4.14) 

(4.15) 
(*.16) 


P  -p  p*c*^i)%;  f- 

0  0 

/ Tl  j2(,)exp(-  -J Tl  p(w»i  -  ’<•'1)5^7  T* 

0  .0  1 

P  ^  -  x)exp(-  1)31  =/Tl  »,(WT1-.,1)5felT  31 
0  0  1 


Then,  using  equations  (4.9),  (4.10),  and  (4.13)-(4.l6), 
we  can  write  the  emergent  intensities  I*(0,+  n)  and  I*(t1,-  |x) 
in  the  forms 
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(4.17)  I*(0,+  -  »/  1  p(m.;t,Ti)^i| 

0 

rT  l  ,  v  B2(t)  d*r 

+  bj  p(H|T1  -  ^^1)w(;iVIiTT  T* 

0 

(4.18)  I*^,-  H)  -  a  f  1  p*(HjTtt1)^|j| 

0 

rTi  */  \  b2(t)  <it 

+  bj  p  (H)Tl  -  ^VwTrf=TT  T* 

0 

From  equations  (4.17)  and  (4.18)  the  emergent  Intensities 
are  given  In  terms  of  the  diffusion  matrix  as  follows t 


/ 1(0*+  M)  \  /Lo(Th 
(4-19)  ,  J-^L,  J 

Vl(Tr-  M-K  * 


where 


(4.20)  10(T)  -  SS^l,  L^t)  -  ^<Tj.  t )  • 

In  equation  (4.19),  the  diffusion  matrix  is  written  in 

the  form 


(4 


•*»  *  -  $  5)  • 

li.  M. 


where  the  reflectance  and  the  transmittance  operators  are 

(4.22)  7^lf(t))  =  f^1  p(^;t,t1)f(t)^-, 

0 

(4.23)  J^f(t))  »/Tl  P#(^t,Tl)f(t)^, 
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(4.24)  ^(f(t))  -/Tl  p'fllJTj  -  ' 

0 

(4.25)  «  /  Tl  p(HiTx  - 

0 

Thus  the  integral  operational  matrix  /X.  given  by  the 

M* 

probabilistic  method  in  the  preceding  paper  [37)  is  now 
derived  from  the  auxiliary  equations  and  the  first  integral 
equation  of  Milne. 

As  a  special  cane,  we  shall  consider  the  diffuse 
reflection  and  transmission  of  a  parallel  beam  of  radiation 
by  a  finite  inhomogeneous  atmosphere. 

Let  a  parallel  beam  of  radiation  of  intrinsic  flux  Fq 
fall  on  the  surface  t  «■  0  in  the  direction  —  Hq,  and  let  a 
parallel  beam  of  radiation  of  Intrinsic  flux  F^  fall  on  the 
surface  t  «=  in  the  direction  +  n^.  Then,  recalling 
equation  (4.20), 

L0(T)  "  TT  exp("  Ll(T)  0  exP(~  jTj-)' 

and  using  equation  (4.19),  we  determine  that  the  emergent 
intensities  are  given  by 

-4^*  M-)  -  |  Fx6(n  -  M-1)exp(—  ~) 

+  x/Tl  Plw’.'iW-M 
0  0 

+  1 j1/  1  p(»in1  -  T,T1)exp(-  j^)^, 

0  1 


(4.26) 
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(4.27)  I(t1#-  d)  -  £  F0«(d  -  dQ)exp(-  Jji) 

+  ^  /  1  p#(^i‘r,r1)exp(-  rr-)^ 

0  ° 

+  ir/Tl  P*(Wt1  ~  T>Ti)exp(- 
o  1 

Writing 

(4.28)  J  1  p(HiT,T1)exp(-  ~-)dT  -  S(t1;hq1h)> 

0  0 

(4.29)  f  1  p(HJT1  -  TiT1)exp(-  jf-)dT  »  T(t1jH1/,h)> 

0  1 

(4.30)  J  1  p*(HjT1  -  T#T1)exp(-  »  S*(t1;h1,h), 

0  1 

(4.31)  J  1  P*(^iT#T]L)exp(-  ~-)dx  -  T*(t1;h0#h), 

0  0 

and  allowing  for  the  principle  of  reciprocity  (cf.  [34]),  i.e», 

-  S(t1;h,^0)  «  S(T1iH0^M.)# 

(4.32)  S*(t1;M.JM-1)  -  S*(T1jH1,n), 

T(T1;n,dQ)  «  T#(T1j^0,d), 

we  obtain 

(4.33)  1(0,  +  d)  -  \  F^d  -  d1)exp(-  ^ ) 

Fn  F,  , 

+  lfil  S(Tl^*^0^  +  ?d  T  (TliM-,dj), 
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T  • 

(404)  H)  -  “  HQ)exp(-  ji). 

P  P 

+  Ijg  T(«r1iH,H0)  +  ^  S*(t1;h,h1). 

The  emergent  intensities  (4.33)  and  (4.34)  are  equal  to 
those  given  in  the  preceding  ijaper  [37]*  Furthermore,  if  the 
optical  properties  of  the  medium  are  constant  throughout  the 
atmosphere,  the  source  functions  p(m-;t,t^)  and  p 
are  equal  respectively  to  p  (M-j^  —  T, t1)  and  p(n;T^  —  t,t^). 
In  this  case,  the  reflectance  and  transmittance  operators  (4.22) 
and  (4.23)  reduce  to  those  yielded  by  Sobolev  [27]. 
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